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Abstract 



We present some observations on the asymptotic behaviour of the coefficients 
in the Laurent series expansion of solutions of the first Painleve equation. For 
the general solution, explicit recursive formulae for the Taylor expansion of the 
tau-function around a zero are given, which are natural extensions of analogous 
formulae for the elliptic sigma function, as given by Weierstrass. Numerical 
results on a solution with pentagonal symmetry are also presented. 



It is known that all solutions of ([T]) are non-classical transcendental functions, in the 
sense that they cannot be written in terms of algebraic or elliptic functions of z, nor 
in terms of classical special functions defined by linear differential equations. The P/ 
equation was discovered in Painleve's classification of second order differential equations 
whose solutions have the property that all movable singularities are poles. For the 
equation ([1]), it is further known that all solutions of this equation are meromorphic, 
with an infinite number of poles in the complex plane. Although Painleve's original 
approach to proving meromorphicity, as outlined in chapter XIV of Ince's book [S], has 
subsequently been found to be flawed, rigorous proofs are now available, with the most 
concise treatment being that of Steinmetz |13j . 

The papers in this special issue are written in honour of Okamoto's insights con- 
cerning the space of initial values for the Painleve equations. The important point 
made by Okamoto is that, in order to consider the solutions of a system with poles, 
such as ([T]), one should enlarge the phase space: as well as regular initial conditions 
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1 Introduction 



The first Painleve equation (P/) is usually written in the canonical form 
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{u{zo) , u' (zq)) e C^, it is necessary to consider the case where u (hence also u') has a 
pole at z = zq] this requires adding points at infinity to the phase space, and leads to 
a sequence of blowups around singular points [2] . 

Boutroux showed that, for suitable z with \z\ large, after the scaling u = z 
and the change of independent variable Z = ^z^^"^, the solution of ([1]) behaves asymp- 
totically like an elliptic function, f/ ~ where p is the Weierstrass function, which 
satisfies the second order differential equation 

P" = Qp' - \92. (2) 

(To be precise, this asymptotics corresponds to g2 = 12 only.) The purpose of this 
article is to present some remarks about the direct comparison between the solutions 
of the Pj equation and the Weierstrass elliptic functions, at the level of exact series 
expansions rather than asymptotics. 

After introducing our conventions of scaling and notation in the next section, in 
section 3 we proceed to compare the coefficients in the Laurent series for the p function, 
which are written in terms of modular forms, with the corresponding expansion of so- 
lutions of ([1]) around a pole. This produces a natural analogy between the lemniscatic 
and equianharmonic elliptic functions, which have the symmetries of the square and 
the hexagon, respectively, and a special solution of Pj with pentagonal symmetry. In 
section 4 we apply the same idea to the tau-function of Pj, and show how both the 
Hirota bilinear form, and an associated equation of degree four, lead to recursive for- 
mulae for the Taylor coefficients in the expansion around a zero. This both generalizes 
the results of Eilbeck and Enolskii on the Weierstrass sigma function [3], and at the 
same time extends some classical expansion formulae due to Weierstrass [12] to the 
case of Pj. As well as being of theoretical interest, we show that these formulae are 
also useful for doing numerical computations. 



2 The equation with parameters 

The programme of classification initiated by Painleve began by considering the neces- 
sary conditions for an equation such as ([1]) to have only poles as movable singularities. 
This requires some local analysis to identify the form that algebraic singularities can 
take at leading order, followed by a power series expansion around any such singular- 
ity to check whether any secular terms can appear; nowadays this procedure is often 
referred to as Painleve analysis (see [6] and references). For the equation ([1]) the only 
algebraic singularities are double poles, ?/ ~ (^ — zo)~^, where zq is the (movable) po- 
sition of the pole. This leading order behaviour then extends to a power series, whose 
resonances (the places where arbitrary coefficients can appear in the expansion) are 
found by substituting 

« ~ (1 + e(. - .o)^) 

into the equation and comparing the leading order terms that are linear in e. In 
this case one finds a quadratic equation in r, with roots —1,6; the root r = — 1 
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corresponds to the fact that the position of the pole at zq is movable, while the value 
r = 6 means that the coefficient of {z — zo)'^ in the expansion should be arbitrary. 
Finally, one must calculate the terms of a full Laurent expansion around the pole up 
to order {z — Zq)'^, and check that no additional counterterms involving log{z — zq) are 
needed for consistency. Hence one obtains a local series representation of the general 
solution, with two free constants corresponding to the resonances. Assuming that such 
a Laurent expansion is convergent in some neighbourhood of zq, this verifies that ([1]) 
has solutions that are locally meromorphic around a pole. Then for Pj, any local 
Laurent expansion around a pole completely fixes a global solution of the equation, by 
analytic continuation. However, purely local analysis is not sufficient to establish that 
all solutions are meromorphic in C, since global estimates are necessary to ensure that 
poles cannot coalesce to form an essential singularity (see e.g. [T3]). 

For what follows, it is worth mentioning at this stage that the details of the Painleve 
analysis for the differential equation (|2]) are almost identical to those for Pj. However, 
if one replaces the 2; in ([1]), or the constant —jg2 in ([2]), by an arbitrary holomorphic 
function f{z), then at the resonance r = 6 one requires f"{zo) = for all Zq, and hence 
/ must be a linear function, in order to have a consistent Laurent series. 

In the rest of the paper we consider solutions of ([I]) with a pole at zq, so it is 
convenient to replace z — )■ z + zq, which adds a constant to the right hand side of 
([1]). Furthermore, for comparison with the equation ([2]) it is useful to rescale u — ?■ 
(— 6A)~^/^n and z — {—6\Y^^z, so that we arrive at the following equation: 

— = Qu-6Xz--. (3) 

The above equation has two parameters, namely A and g2, and any solution of the 
Pj equation in the canonical form ([1]) with a double pole at 2; = 2:0 corresponds to a 
solution of ([3]) with a pole at z = 0, for a suitable choice of the constant g2, with any 
A 7^ 0. Moreover, when A = the equation reduces to (E]). 



2.1 Recursion for Laurent coefficients 

To determine the local Laurent expansion of the solution of ([3]) with a double pole at 
the origin, we insert the expression 

u{z) = J2cnZ'"\ Co = l, (4) 

n=0 

into the equation, from which subsequent coefficients are determined by the recursion 

n— 1 ^ 

(n + l)(n - 6)cn = 6 ^ CjCn-j - -g2SnA - 6A(5„,5, n > 1. (5) 
i=i 

The above relation determines all the coefficients ci, . . . , C5 uniquely, but cq is arbitrary 
(corresponding to the aforementioned resonance). Once ce has been fixed, then all Cn 
for n > 7 are all uniquely determined by ([5]). We find 

ci = C2 = C3 = 0, C4 = 1^, C5 = A, cg = C7 = 0, . . . , (6) 
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where the value of cq is given in terms of an additional parameter (73, which is introduced 
for ease of comparison with the elliptic case. The following result is immediate. 



Lemma 2.1 For n > 1 the coefficients of the expansion Q), with cg = (73/28, are 
polynomials in g2, A and g^, 

Cn = P„,(fi'2, A,5(3), 

where P„ has rational coefficients and is weighted homogeneous of total degree n in 
these arguments with weights 4, 5 and 6 respectively: 

Pn(C'^72,C'A,C'(73) = CPn((72, A,(73) VC G C*. 

The polynomials P^ for n = 1, . . . , 7 can be read off from ([6]). Note that, a priori, 
they must be identically zero for n = 1, 2, 3, 7, since there are no non-zero homogeneous 
polynomials with these weights; but they are non-zero for all other values of n. Here 
we list the next few non-trivial ones: 

- 1200' ^' - 22 + 280 J ' 140' 

1 /(7i on _ngl\ _ 1 (mg2\' , g,gl 



208 V 750 49 y 49000 44 V 500 4200^ 

These polynomials have some interesting properties. As will be explained in section 
3, for A = the polynomials P2n{92,0, ds) are polynomials in g2 and g^, which define 
modular forms: up to overall scaling they give the Eisenstein series of the elliptic curve 

= 4a;3 _ _ 

In fact, viewing each P„ as a polynomial in A with coefficients in Q[g2,g3], each coef- 
ficient is weighted homogeneous in (72,(73 with some weight, so also defines a modular 
form. Thus we refer to the P„ as modular polynomials. 

Henceforth we let u{z; g2, X, gs) denote the meromorphic function defined by the 
series @. Our main task in section 3 will be to consider the asymptotic behaviour of 
the coefficients c„ in this series as n — 00, and the way this depends on the parameters 
(72, A, (73. Before doing this, we introduce some objects needed in section 4. 



2.2 Hamiltonian and tau-function 

Okamoto discovered the Hamiltonian structure of the Painleve equations [10]. The 
equation can be written as a Hamiltonian system 

, dh , dh 



where the Hamiltoniaiil] is 



1 2 3 1 . 1 

h= -V - 2m + -g2U + <6\zu + -5(3. 



^The constant 53 is included for later convenience, but does not affect Hamilton's equations. 
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Taking the total derivative of h with respect to z we have 

h' = ^ = QXu. (9) 
oz 

The tau-function associated with a solution of ([3]) is a function t = r{z) such that 

« = -^log^- (10) 

Observe that the preceding relation only defines t up to gauge transformations 

T —7- exp(a2; + b) t, a, b arbitrary, (11) 

and the fact that u is meromorphic means that t is holomorphic, with t having simple 
zeros wherever u has double poles. 

We also note that from Q, up to the addition of a constant we have 

/i = -6A-^logT. (12) 

dz 

The constant of integration has been fixed so that when A = we have h = 0, and 
the right hand side of ([8]) corresponds to the first order differential equation for the 
Weierstrass p function, i.e. 

{p'f = Ap^-g2p-g3, (13) 

for u = p and v = p'. This choice of integration constant will also be convenient for 
the analysis of the tau-function when A ^ 0. 



3 Asymptotics of coefficients in the Laurent series 

We now describe the behaviour of the coefficients c„, which satisfy the recursion 

CkCn-k, n>7 (14) 



b 



{n + l){n-Q) 



k=l 



for a suitable choice of initial values. To begin with we consider the classical case A = 0, 
with u{z; g2, 0, g^) = p{z; g2, gs), for which ci = C2 = C3 = C5 = in (111]). 

3.1 The elliptic case 

The Weierstrass p function can be defined as the unique solution of the first order 
differential equation ( IT3l) having a double pole at 2; = 0. However, in all standard 
treatments of elliptic function theory (see e.g. chapter XX in [E]), it is usually defined 
by its Mittag-Leffler expansion, that is 

p(z;g2,g3) = ^+ Yl ( (z-ny --^) ' 



!^6A\{0} 
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where A denotes the period lattice, consisting of all periods of the form = 2mooi+2nuj2 
for {m,n) G Z^, generated by the two independent periods 2uji, 2uj2- The uniform con- 
vergence of the expansion (fT5|) is guaranteed by the fact that the power sums ^^^.^q 
are absolutely convergent for all /i > 2, and expanding around z = gives the series 
(jl]) with Co = 1, Ci = C2 = and 

c„ = (n- 1)^1]-", n>3. (16) 
Thus it follows from f lTSj) that p satisfies the differential equation f|T3|) with 



Conversely, given a solution of f[T3|) with parameters g2,g3, two independent periods 
are obtained from the integrals 2uj = j^ ^, j = 1,2 around cycles Ci, C2 that generate 

the homology of the elliptic curve ([7]). (This is in the generic case — 27 ^ 0; 
otherwise (fT3l) is solved in elementary functions.) 

Observe that, from the symmetry of the lattice A = —A, where the poles of p{z) 
are situated, the odd index coefficients C2n+i = for all n > 1 (hence p is an even 
function). The even index sums (fTB]l can be rewritten as 

2n — 1 

= ^2uj^fn '^2nir), (17) 

where the function G2n{'T) is the Eisenstein series, a modular form of weight 2n corre- 
sponding to a normalized lattice with periods 1, r, with t = 

G2n{r) = Yl 



in + ToY"^ 

(p,g)eZ2\{(0,0)} 

The periods 2cji, 2u}2 can always be ordered so that Imr > 0, and by choosing 
|2a;i| to be minimal we also have |r| > 1. Then the series dl]) has radius of convergence 
|2a;i|, and the asymptotic behaviour of C2„ depends on r. There are two possibilities: 



|r| > 1: lim„^ooG'2n(r) = 2; 

|r| = 1: G2n{T) has no limit as n — )■ oo. 



Note that we can always choose r so that it lies in the fundamental domain J-" = 
{|t| > l}n{|Rer| < \} ([12], Proposition 1.5). In the first situation, ±1 are the 
only normalized periods that lie on the unit circle, so all of the terms in the sum 
(|T8l) vanish as n — )■ oo apart from the contribution from (p, g) = (±1,0). For the 
second possibility, r = e'^ lies on the lower boundary of J-", and (by using the modular 
transformation r — )■ — 1/t) it can always be assumed that vr/S < 6 < 7r/2; then 
9 = n/3, with six normalized periods on the unit circle (the sixth roots of unity), or 
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/ i 




1 


5.86303169342540159797 


2 


6.00963997169768048102 


o 

o 


5.9997183563705/593409 


4 


6.00001 164 <^0<^9<^ <^9M4o5 


r 



5.99999950 /^4o55330152o 


u 
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11 


5.99999999999999892076 


12 


6.00000000000000003997 


13 


5.99999999999999999851 


14 


6.00000000000000000005 



Table 1: The oscillating values of G^ni^ ^) • 



otherwise there are four periods ±1, ±r with modulus one. The extreme cases are the 
so called equianharmonic case [9 = 7r/3), and the lemniscatic case {6 = 7r/2). 

Equianharmonic case: This corresponds to (72 = in (fT3|l . Without loss of generality 
(by rescaling) one may set gs = 1 and find half-periods 

= — a;2 = e3a;i, (19) 
y ^ V 4x3 - 1 4:71 

which gives r = eT. Then from the recurrence (1141) with C4 = it follows that 
G2n(e^) = unless n = mod 3. The Laurent series for the p function becomes 

Pfe0.1) = l + f ^g^Ga^e*).""- (20) 

n=l 

with uji as in f lT9|) . There are six normalized periods , j = 0, ... ,5 on the unit 
circle, so in the limit of the Eisenstein series flTSj) . only six terms survive to yield 

lim Gdnie^) = 6. 

Some values of these modular functions are given to twenty decimal places in Table [H 

Lemniscatic case: This case corresponds to = 0, and (upon scaling so that 
g2 = 4) the half-periods are found as 

dt I , , 
uj, = -j=== = -5(1/4, 1/2), U2 = ia;i, 

whence r = i. The fact that cg = in (IHj) now implies that G'2n(i) = unless n is 
even, so that the series (jl]) in this case takes the form 

1 °° 4„ _ 1 

K-;4,0) = - + 5^-^G4„(i)z--. (21) 

n=l ' 
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3.15121200215389753821 
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4.25577303536518951844 
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4 


A ni no ono [X no /I o OT 
4.UiooyoUooUzoUz4oooo ( 





o nn/^nn/^ Tlx o 1 '7^ oon [X txn [X '7 

o.yyouyo /^ooi / ozoyooyo ( 
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11 


3.99999904632591103400 


12 


4.00000023841859318284 


13 


3.99999994039535611558 


14 


4.00000001490116124950 



Table 2: The oscillating values of G4n(i) • 



The coefficients can also be expressed in terms of the Hurwitz numbers [H [3 E] 
(elliptic analogues of the Bernoulli numbers) which are given by 

^" = (22) 
From (fT4l) also follows the recurrence found by Hurwitz: 



As can be seen from Table [2l the limit of the Eisenstein series ( !T8|) in this case is 



(23) 



lim G4n(i) 

n—^oo 



3.2 The general case A 7^ 

All solutions of Pj are known to have order of growth 5/2 [T3], and this is the same as 
the infimum of the values of /i such that the power sums over non-zero poles, Xln^^o ^~^) 
are absolutely convergent (|I5], chapter VIII). As a consequence, the solutions of ([3]) 
with A 7^ admit a Mittag-Leffier expansion of the same form f lTSj) as for p{z): 

^' = i + E ( - ^) (24) 

poles 

(cf. [13], Theorem 7.1). The main difference with the elliptic case is that the poles of 
u no longer lie on a lattice, and suitable analogues of elliptic integrals to determine the 
positions of the poles are unavailable. 
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Nevertheless, the formula (flGjl for the coefficients in (j4]) still holds in general, and 
as a consequence, from (Q, C3 = = cy gives Xln^o^"^ = = Xlri^o^"'^' "while 
generically one expects that c„ 7^ for n > 4, n 7^ 7 (and similarly for the corresponding 
power sums). Following the elliptic case, by picking a non-zero pole fi* such that 
is minimal, one may write 

* n^o 

where the sum is over the non-zero poles Cl = Q/Q^, of the rescaled solution 

uiz; 5(2, A, ^3) with g2 = ^t92, A = fi^A, ^3 = illg^. 

Generically, for the rescaled solution, there should be only one pole with modulus 
1 (corresponding to Q^), with all other non-zero poles lying outside the unit circle. 
Assuming this to be the case, it follows that 

lim F„ = 1, whence = lim — —. 



n—^oo 



Thus, in the generic situation, we find that numerical iteration of ( 1T4|) allows the pole 
nearest to 2 = to be calculated very efficiently. 

The question of precisely when a non-generic solution can arise, with two or more 
poles Q having the same modulus seems to be a very difficult one. However, there 
is one situation where we know this to be so, namely when (72 = = ^'3: this is the 
pentagonal case. 

Pentagonal case: When g2 and gs both vanish, it follows from (|5]) that c„ = unless 
n = mod 5, so that the Laurent series (j4]) for u{z) = u{z; 0, A, 0) takes the form 

u{z) = -2+J2^-^F^X-z'-^ (26) 

^ n=l ^ 

where we set 

(5n-l) - _ 

Above we introduced the coefficients f„ which satisfy the recurrence 

6 "^"^ 

y^^fcfn-fc, n>2, with vi = l. (27) 



(5r2 + l)(5n-6) 



k=l 



In this case, the non-zero poles of the solution ( 1261) lie on regular pentagons centred 



at 2; = 0: when is a pole of u, then so is e f2 for j = 1, 2, 3, 4. We have F„ = if 
n 7^ mod 5, and from (1251) . with 5 poles on the unit circle, 

lim = 5. 

n— >oo 
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4.99999999957591996469 


12 


5.00000000005151463070 


13 


4.99999999999374379484 


14 


5.00000000000075986460 



Table 3: The oscillating values of Fy 



At n = 30 we find a value of 7 apparently correct to 23 decimal places: 

7 = lim / " ^ 18.32138268472483887119960.... (28) 

n-5.oo [pn — l)Vn+l 

Given this accuracy in 7, the quantities Fn = F^n can then be calculated, and should 
converge exponentially fast to the value 5; see Table |3l The value 7^/^ ~ 1.788923 gives 
the non-zero real pole of the solution u{z; 0, 1, 0) closest to the origin, and the radius 
of convergence of the series (!26|) : for an independent verification of this see section 4. 

We conclude this section with some observations on the cases (72 = and = 0. 
3.2.1 The case 5-2 = 

To consider the form of the quantities c„ when 52 = 0, it is convenient to use the 
parameter a = 53/28. Thus we have the recursion (fT^ with cq = 1, Ci = C2 = C3 = 
C4 = 0, C5 = A and cq = a, and find that the structure of the iterates depends on the 
index mod 5, so that 



J2 4ni, a5-+PA"-6"^-^ p = 0, . . . , 4, (29) 

m=0 

for certain rational numbers Cg^^^, which satisfy quadratic recurrence relations. 

To begin with, we examine cfj, which are the coefficients of a^. Upon setting 
Vn = cfj, we see that Vn satisfies fl27|) . The next terms to consider are cfn+i^ i-^- the 
coefficients of a^, and upon setting Wn = cfn+i 



12 



Wr, 



(5n + 2)(5n-5) 



''^^WkVn-k, n>2, with Wi = l. (30) 



10 



Observe that, supposing the sequence of Vn to be known, the above recurrence is linear 
in the unknowns (unhke the recurrence (1271) for the Vn themselves). This means that 
the generating function for Wn satisfies a second order hnear ODE, within which the 
generating function of the u„ appears as a coefficient. 

As we shall now see, this property extends to the recursive generation of all the 
sequences c^].p when {m,p) ^ (0, 0). The linear recurrence solved by the general term 

4n+p be written explicitly as 



n— 1 / m p 1 4— p 

■^'f I •"If 

k=i \j=o e=o j=o 1=1 



A:=l \i=0 

where A'n.p is given by 



^-^v=T^-——^rr^ Try P = 0,1, 2, 3, 4. (32) 

(on + p + 1 j (on + p — 0) 

The initial conditions for the equations (I3T]) are that 4""^ = for 1 < A; < 9 and for all 
m > 0, except for Cg''^ = Cg°'* = 1. From (PT]) it appears that all the nonlinearity in the 
problem of determining these coefficients is moved into finding the solution of (12 7p . 

For the recurrence (^7^ , note that it is sufficient to determine a single no n- vanishing 
solution, since all other solutions can be obtained by the rescaling Vn A"'Vn, with 
A arbitrary. The Vn correspond precisely to the coefficients in the Laurent expansion 
(l26l) of the solution u{z) with pentagonal symmetry that was discussed previously. To 
see an example of how this is related to the other coefficients, we consider Wn = cfn+i 
once more, and introduce the generating function G{x) = Yl'^=i WnX'^~^. From f pOj) it 
follows that G satisfies the ODE 

^ (fa ^ 12 (IG _ 12 
dx"^ 5 dx 25 

where tplx) = Xl^i 'i^nX"'~^- In fact, up to some scaling and shifting, iIj{x) is just given 
by the pentagonal solution u; to be precise, 

x\^ ( fx 



X 



In a similar way, via ( l3T|) . the generating functions for the other terms c"^\p are related 
to each other: for example, that of €^^^2 is related to the generating functions of the 
sequences cf^ and c^n_^_i'., and that of Cg^^ is related to the generating functions of the 
sequences cf^^p, p = 0, . . . , 4, and so on. 

3.2.2 Case 5-3 = 

For the case 53 = 0, it is helpful to introduce the parameter /3 = 52/20, and take the 
solution of IHM with cq = 1, Ci = C2 = C3 = Ce = 0, C5 = A and C4 = /3. The coefficients 
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of (jlj) now have the structure 



m=0 



for some rational numbers Cg^lp. 



J2 ct-p /3'^+^X^-'^-P, p = 0, . . . , 4, (33) 



The quantities CgJ, which appear as coefficients of are precisely the same as 
the numbers Vn = cf^ found previously, since setting /3 = just gives the pentagonal 
solution u with the expansion fj26|) . The next terms to consider are w„ = c^^n-i^ ^^e 
coefficients of which satisfy the recurrence 



12 



Wr. 



(5n)(5n-7) 



^^WfcWn-fc > 2, with wi = 1. (34) 



Just as for the case g2 = 0, once Vn is given, this recurrence is linear in the unknowns 
Wk- Analogously to equation (13T|) . it is possible to write down the general recurrence 
solved by the term Cg^lp, the result being 

/ m p n—1 m—1 4— p n \ 

— r^n-p I ^ ^ ^ Hk-i^b{n-~k)+i~p ^ ^5k-i-p^^{n-k)+l I ' l"^^^ 



0=0 e.=Q k=l j=0 l=\ k=l 



■'bn—p 



where again Kn,p is given by ( 1321) . but above it appears with p — )■ —p compared with 
fl3Tl) . The formula fl35|) holds for n > 2, and the initial conditions are given by c^"*'' = 
for 1 < < 6 and for all m > 0, except for £4°^ = £5°^ = 1. Similarly to the situation 
for (72 = 0, linear ODEs for the generating functions of the rational numbers Cg^l^ for 
(m,p) 7^ (0,0) can be constructed recursively, once the numbers Vn are known. 



4 Expansion of the tau-function 

The tau-function for ([3]), related to u by ffTOl) . satisfies a fourth order differential equa- 
tion which is homogeneous of degree two. It is written in Hirota bilinear form as 

D^T-T- (12A^ + ^2)t2 = 0, (36) 

where denotes the Hirota derivative defined by 

Painleve analysis can be applied directly to the equation (136|) . expanding around a 
simple zero at 2; = 0, corresponding to a double pole in u there. Seeking resonances by 
taking t ~ 2; + ez^ yields r = —1, 0, 1, 6: the value —1 is the movable position of the 
singularity, as usual, and the values r = 0, 1 correspond to the two free parameters a, h 
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for the gauge transformations ffTTj) . which leave the equation (136!) invariant; this leaves 
only r = 6, which is equivalent to the freedom to choose cq (or the parameter g^) in 
the Laurent series (jlj). 

We would like to determine the Taylor series of the tau-function around z = 0, with 
t(0) = 0, since (due to the fact that t is an entire function), this provides a global 
representation of the solution of ([3]), via the formula ( fTOi) . If we fix the gauge, we can 
always choose the coefficient of 2 to be 1, and set the coefficient of z"^ to be 0, which 
results in a series of the form 

oo 

T(z) = z + 5^a^"+^ (37) 

n=2 

This is the analogue of the power series for the Weierstrass sigma function (t{z), which 
satisfies the bilinear equation fl36l) with A = 0. Using the same method as in [3], with 

Dtz^ . z' = b,,,z=^'-\ where b,,, = 4! g(-l)^ ( ( 4 " ' 
it is straightforward to establish the following. 



Theorem 4.1 The coefficients in the expansion ( [^f7| j of the tau-function belong to 
727^,(73], being uniquely determined by the recursion 



n— 1 ^ n— 4 n— 5 

n{n'^ -l){n- 6)C„ = "2 &j+i,n-i+iCjC'n-i + -92 ^ Cfin-i-j + 6A ^ CjC„„5_j 

j=l j=0 i=o 

subject to fixing Co = l,Ci = 0,Cq = —gs/SAO. Each Cn is a weighted homogeneous 
polynomial of total degree n in the arguments (72, A, (73 with weights 4, 5, 6 respectively. 

Remark 4.2 The choice of gauge is specific to the zero at z = 0. If we let t:{z; 92, A, g^) 
denote the function given by (31)' ^'^^ then expand around another zero at z = VL ^ {), 
we obtain the formula 

x{z- g2, A, gs) = Ae'"^'-'''>x{z - Q; g2 + 12\n, A, ^3), 

where in principle g^, A = T'(f2) and B = ^r" (^l) / x' (il) all depend on g2, A, (73 and 
The quasiperiodicity of a{z) under shifting by a period is a special case of this. 

For the case of the sigma function in [3] it is noted that the value of the coefficient 
Cq must be given appropriately in terms of gs in order to be consistent with the Laurent 
series for p satisfying (fT3|) . The same holds for A 7^ 0, upon requiring consistency of 
(1371) with (jl]), as is seen by noting that the expression ([8]) for the Hamiltonian gives a 
first integral for the bilinear equation (!36|) . taking v = u' with u given by (fTO!) and h 
given by fll2p . Then the formula ([H]) can be rewritten as a third order equation for t 
that is homogeneous of degree four, namely 

t2(t'")' - 6tt'tV" + 4(t') V" + 4t(t")=' - 3(t't")' 

-g2X^{xx" - (t')2) - l2\[z{xH" - ixx'f) - t^t') + g-^x^ = 0. ^^^^ 
This equation immediately yields another recurrence for the coefficients of t. 
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Figure 1: The poles of the solution of ([3]) with 32 = ffs = 0. 



Theorem 4.3 Subject to fixing Cq = 1,Ci = 0, the coefficients in the expansion \31\) 
are uniquely determined by the recursion 

(4) (5) (6) 

n{n^-l)Cn = - ^ bj^k,i,m CjCkCjCm + 92y^, +12A -g^ 

j+k+l+m=n 
l<j,k,£,m,<n—l 

where 

hk/,nr = j{j + l)ik + l)({j-l){k{k-6i-7)+4{i+l){m + l)) + k{i+l)i4i-3m-3)y 
and 

(4) 

J2= E U + l){j-k-l)C,C,C,Cr., 

i+k+l+m=n—A 
0<j,k,£,m<n—A 

(5) (6) 

j+k+£+m=n—5 j+k+l+m=n—Q 
0<j,k,£,m,<n—5 0<j,k,£,'m<n—6 

Both of the recursions for C„ given above can be iterated very rapidly to produce 
polynomial approximations to the tau-function, for any choice of the values g2,^,g3- 
As an example, to generate Figure [T] we took the values g2 = ^ = 93 and A = 1, 
corresponding to the solution with pentagonal symmetry that was considered in the 
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previous section, and calculated the first 101 non-zero terms in the expansion (j3 
which contains only terms of the form 2^""*"^, beginning with 

Iz^^ z^^ 83^21 

~ ^ ~ 20 ~ 26400 ^ 1232000 ^ 117976320000 ^ "" 

In the figure we have plotted the roots of the truncated series, which is of the form 
-P50i(^) = ^ Aoo('?^)! where Pioo is a polynomial of degree 100 m. w = z^ . The roots 
of the truncated series approximate the zeros of T with increasing accuracy as more 
terms are added. The smallest root of Pioo is a real number 7, agreeing with the digits 
of the numerical value (!28|) to the same accuracy as before. In fact, it appears easier 
to get improved approximations to this value by truncating the Taylor series than by 
calculating the ratios of the coefficients c„ in the Laurent expansion. 

Before concluding this section, we present another representation for the series 
(l37Il . which displays some interesting arithmetical features of the coefficients. Note 
that, since the C„ are weighted homogeneous, the tau- function can be written in the 
form of a triple sum 

1 4£+5m+6n+l 

E A...(;..)'(6Ar(2,3r (39) 

^,m,n>0 ^ ' 

for certain rational numbers Ai^m.n- The above formula is motivated by the result of 
Weierstrass [I6], who showed that the sigma function can be expressed as 

1 ^4m+6n+l 

a{z)= a^,A-g2)"'{2g3rj-—-——. 
^-^ 2 4m + 6n + 1 ! 

m,n>0 

From comparison of the two series it is clear that upon setting A = we have 

am,n = ^m,0,n Vm, K > 0. 



Weierstrass used ([2]) and f[T3|) . together with various modular relations, to show that 
the sigma function satisfies two linear partial differential equations, from which he 
obtained a linear recurrence relation for the coefficients a^.n- The first of these PDEs 
just follows from Euler's theorem on homogeneous functions, and thus extends to the 
tau-function when A 7^ also: 

d , d d d \ 

4^27^ + 5A— + 6^3^ 2— + 1 T = 0. 

092 oX dgs dz J 

(Note that z has weight —1, so overall t has the same weight.) However, the other 
PDE that Weierstrass found for cr(z) is of second order in z (containing also first 
derivatives with respect to 92,93)', it is equivalent to the fact that the corresponding 
elliptic theta function satisfies the heat equation. For the tau-function of Pi we do not 
expect any linear PDE of this kind. Nevertheless, the Hirota bilinear equation fl36|) 
provides recursive relations for the coefficients Ai ^ n- 
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Theorem 4.4 Subject to fixing Ao,o,o = 1? ^0,0,1 
completely determined by the recursion 



—3, the coefficients in (E 



are 



s{s-l){s-2){s-l) 



E' 



1 

'2 ^ Si!S2! 

mi+m2=m 
n\+n2=n 



+ E 

h+l2=i 
mi+m2=m—l 
n\+n2=n 



,mi ^^2 1^2 '^2 
Sl\S2\ 



E Al-^ .mi.ni ,m,2 ,112 

Sl!s2! ' 

ei+e2=i-i 

m\+m2=m 
ni+n2=n 



where s = M + 5m + 6n + 1, = Mj + Sm^ + Grij + 1 for j = 1, 2, and denotes 
that the terms with {£i,mi,ni) = (0,0,0) or {£,m,n) are omitted from the sum. 

Remark 4.5 Another recursion for A^^rn^n can be obtained from / fS'Sj) . 

The first few coefficients are given below in the form of 3 x 3 matrices M*^™^ whose 
(j, k) entry is Aj^i^rn,k-i for m = 0, 1, 2: 



1 -3 
-1 -18 



-54 
4968 



-9 513 257580 
M(2) = 



-294 
18774 



144144 
15053040 



1112436 3160803600 



-6 -216 89424 
-84 18720 5786640 
1650 1358640 1168920720 

47585880 
22914336240 
-2734614623160 



Due to recent results of Onishi [TT], it is known that the coefficients a™ « in the 



expansion of the sigma function are all integers, and analogous results have been proved 
for sigma functions of some higher genus curves [1] . This suggests 



Conjecture 4.6 The series ^3^) has Ae^m,n G Z Wi,m,n > 0. 
We have verified this conjecture for the first few hundred coefficients. 



5 Conclusions 

The properties of the exact series expansions for the solution of Pj and its tau-function 
are tantalizingly close to those of the analogous Weierstrass functions. As well as their 
potential uses for numerical calculations, it would be interesting to use the recursions 
in Theorems 14.11 and 14.31 to prove directly that T is holomorphic (hence providing yet 
another proof of the fact that all solutions of Pj are meromorphic). Finally, we should 
mention that all of the Painleve equations give deformations of elliptic functions, and 
they all have associated tau- functions and bilinear equations [5] , so the same methods 
could be applied to Pu^vi- 
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